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Abstract. We present the results of a detailed timing 
analysis of observations of Cen X-3 taken by the University 
of Durham Mark-6 imaging atmospheric Cherenkov tele¬ 
scope in 1997-1999. The presence of a TeV 7-ray signal at 
the overall > 4.5cr significance level in the ‘fully cut’ image 
selected data, as reported earlier, is confirmed. A search 
for possible modulations of 7-rays with the pulsar spin 
period Pq ~ 4.8 s was performed by the step-by-step ap¬ 
plication of image parameter cuts of gradually increasing 
hardness. The data of each of 23 days of observations have 
not revealed any statistically significant Rayleigh power 
peak, except for 1 day when a peak with a chance prob¬ 
ability p = 6.8 X 10“^, was found in ‘soft-cut’ data sets. 
This modulation, if real, is blue shifted by 6.6 msec (> 10^ 
kms“^) from the nominal second harmonic of the X-ray 
pulsar. Taking the large number of frequency trials into 
account, the estimated final probability of such a peak by 
chance still remains < 10“^. Bayesian statistical analysis 
also indicates the presence of such modulations. We show 
that the behaviour of the Rayleigh peak disappearing in 
the fully cut data set is actually quite consistent with the 
hypothesis of a 7-ray origin of that peak. No modulation 
of the VHE 7-ray signal with the pulsar orbital phase is 
found. 

In the second part of the paper we consider differ¬ 
ent theoretical models that could self-consistently ex¬ 
plain the existing data from Cen X-3 in high-energy (HE, 
E > 100 MeV) and very high energy (VHE, E > 100 GeV) 
7-rays. We propose on the basis of the energetics re¬ 
quired that all reasonable options for the 7-ray produc¬ 
tion in Cen X-3 must be connected to jets emerging 
from the inner accretion disc around the neutron star. 
One of the principal options is a large-scale source, with 
Rs ~ 10^^ — 10^^ cm; this assumes effective acceleration 
of electrons up to ~ 10 TeV by shocks produced by in¬ 
teraction of these jets with the dense atmosphere of the 
binary. It is shown that such a quasi-stationary model 
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could explain the bulk of the 7-radiation features observed 
except for the 7-ray modulations with the pulsar spin. 
These modulations, if genuine, would require an alterna¬ 
tive source with Rs 10^^ cm. We consider two principal 
models, hadronic and leptonic, for the formation of such 
a compact source in the jet. Both models predict that the 
episodes of pulsed 7-ray emission may be rather rare, with 
a typical duration not exceeding a few hours, and that gen¬ 
erally the frequency of pulsations should be significantly 
shifted from the nominal frequency of the X-ray pulsar. 
The opportunities to distinguish between different models 
by means of future 7-ray observations of this X-ray binary 
are also discussed. 

Key words: acceleration of particles - radiation mecha¬ 
nisms: non-thermal - X-ray binaries: individual: Cen X-3 
- gamma-rays: theory - gamma-rays: observations 


1. Introduction 


Centaurus X-3 has been one of the prominent galactic 
sources of hard radiation since its discovery as one of the 
first cosmic X-ray sources (Chodil et al. 1967 ), and was 
the first X-ray pulsar to be discovered in a binary system 
(Giacconi et al 1971, Schreier et al. |1972 ). All the basic 
parameters of this archetypal high mass X-ray binary are 
well known. The pulsar has a spin period Pq « 4.8 s and 
an orbital period Porb ~ 2.1 d, with a gradual shorten¬ 
ing (i.e. ‘spinning-up’) of both periods in time, and with 
a deep eclipse of the X-ray source at the orbital phases 
—0.12 < (j) < 0.12 (for reviews see Joss & Rappaport 


1984, Nagase 1989). The X-ray luminosity of the pulsar is 
very large, reaching Lx ^ 10^®ergs“^ in the ‘high’ state 


(e.g. White et al. 1983| , Burderi et al. 2000), which im¬ 
plies a massive accretion of material onto the neutron 
star from the optical companion (Pringle & Rees 1972| , 
Lamb et al. 19731 ). The optical companion (V779 Cen) 
was discovered by Krzeminski ( |1974 ) and has been iden¬ 
tified as an evolved 0-type star with surface temperature 
T > 3 X 10'^ K at a distance ^ 8 kpc from the Sun, and 
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a bolometric magnitude of Mboi ~ —9 (Hutchings et al. 


1979). The masses of the stars in this binary are estimated 


as Mn ~ 1.2 Mq and Mo-star — 20 Mq. The value of Porb 
suggests the separation between the centres of the stars 
to be a = 1.3 x 10^^ cm and the radius of the massive star 
filling its Roche lobe Rq = 8.6 x 10^^ cm (see e.g. Clark 
et al. 1988, Ash et al. 19991 ). 


Cen X-3 is also known as one of those X-ray bi¬ 
naries from which y-ray signals have been reported. In 
the domain of high energy (HE) y-rays, conventionally 
E > 100 MeV, the y-ray flux during 2 weeks of obser¬ 
vations of Cen X-3 in October 1994 by EGRET was at 
the level F{> 100 MeV) = (9.2 ± 2.3) x 10“^phcm“^s“^ 
(Vestrand et al. 1997 ). This excess flux was significant at 
a 5(7 level, which is generally considered as a reliable de¬ 
tection with the EGRET instrument. The timing analysis 
has shown a significant modulation of the signal with the 
pulsar spin, precisely at the contemporaneous frequency 
of X-ray pulsations measured by BATSE, with a chance 
probability estimate based on the Rayleigh test statistic 
of about 1.6 X 10“^. No modulation with the orbital phase 
of the pulsar was found: the y-ray signal seems to be quite 
homogeneously distributed throughout the entire orbit, 
with 68 of the 264 HE y-rays from Cen X-3 being de¬ 
tected in the pulsar eclipse orbital phase \(j)\ < 0.12. This 
suggests a production site of the y-rays far away from the 
pulsar. The y-ray signal was not found in the data of other 
observing periods of Cen X-3 by EGRET, which suggests 
a significant variability of the HE y-ray source on a time 
scale of several months (see Vestrand et al. 1997 for de¬ 
tails). 

In the domain of very high energies (VHE), con¬ 
ventionally E > 100 GeV, sporadic y-ray signals from 
Cen X-3 had been earlie r repo rted by the University of 
Durham (Brazier e t al. 1990a ) and the Potchefstroom 
groups (North et al. 1990 ). Note that the early detections 
of VHE y-rays from X-ray binaries in the 80’s had been 
carried out by non-imaging Cherenkov telescopes which 
had rather poor sensitivity. The y-ray signals were there¬ 
fore extracted mostly on the basis of the timing analyses, 
and are subject to controversy concerning their reliabil¬ 
ity (see Weekes 19921) . It is worth noting in this regard 
that Cen X-3 is the first and until now the only X-ray bi¬ 
nary to be detected as a source of if > 400 GeV y-rays by 
a contemporary imaging instrument, the Durham Mark- 
6 telescope, as we have reported earlier (Chadwick et al. 
199S, ^000 ). The excess (y-ray) signal has been found in 
the ‘ON-source’ data during each of the 3 years of ob¬ 
servations from 1997 to 1999. The estimated mean y-ray 
flux is about E{> 400GeV) ~ 2.8 x 10“^^ cm“^s“^, at 
the significance level for the entire data set of 4.7 (t. No 
significant modulations in the combined ON-source data 
with either the 4.8 s pulsar period or 2.1 d orbital period 
of the binary were found, but these data were analyzed 
after application of the image cut procedures (Chadwick 


et al. 200C) which we shall see later can be very counter 
productive in a periodicity search. 

In this paper we present the results of a more detailed 
timing analysis of observations of Cen X-3, which in par¬ 
ticular includes a step-by-step search for possible short¬ 
term episodes of y-ray emission modulated with the spin 
of the pulsar, as well as a search for a modulation of the y- 
ray signal with the pulsar orbital phase (Section 2). Then 
we carry out in Section 3 a detailed theoretical study of 
the consistency of different models for the production of 
y-rays in Cen X-3 with the experimental data currently 
available in both HE and VHE domains. Finally, in Sec¬ 
tion 4 we summarize the results of our current study, and 
discuss some possible tests for future y-ray observations 
that could help to distinguish between different models of 
y-ray production in this X-ray binary. 

2. Timing analysis 


As we have reported earlier (Chadwick et al. 1998, 2000), 
the entire set of data of observations of Cen X-3 by the 
University of Durham Mark 6 telescope (see Armstrong 
et al. 1999 for a description) from Narrabri, Australia, 
consist of pairs of 15 minute segments of alternately ON- 
and OFF-source observations. A total of 129 ON-source 
segments have passed the relevant consistency criteria in 
order to be accepted for further analysis, yielding about 
32hrs worth of data. The data have been accumulated 
during 23 days of observations over 3 years from 1997 to 
1999. After application of a set of shower image parame¬ 
ter cuts (see Chadwick et al. 2000 ), an excess of events in 
the on-source data, indicating a possible presence of y-ray 
signal, is found in the data from each individual year. In 
Table 1 we present the results of the analysis of the data 
for all of the 23 days of observations. Note that the signif¬ 
icance 4.34(7 in the overall data set for the image bright¬ 
ness parameter brightness > 800 digital counts, which 
corresponds to the Monte-Carlo estimated y-ray energies 
E > 400 GeV, is somewhat lower than 4.7 (t as reported 


previously (Chadwick et al.200C) because two days worth 
of data were subsequently corrected for a problem with 
one of the imaging PMT’s. 

For the periodicity analysis the air shower arrival times 
were first corrected to the reference frame of the solar 
system barycentre in order to exclude the timing effects 
connected with the motion of the Earth. 

2.1. Orbital phase modulations 

Important information about a plausible site and mecha¬ 
nism of gamma-ray production in X-ray binaries can be 
derived from studies of possible correlations of the y-ray 
fluxes with the orbital phase of the pulsar. This can be 
particularly informative in the case of Cen X-3 because it 
is an eclipsing binary, and because its optical c ompan ion 
is a very luminous star. As shown by Bednarek (200C|), in 
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Orbital Phase 



Orbital Phase 

Fig. 1. ON-OFF excess, in terms of ‘ 7 -rays per cosmic 
ray’, in the fully cut data sets with 2 different brightness 
parameters, plotted with respect to orbital phases of ob¬ 
servations in 1997-1999. 


the case of 7 -ray production near the orbit of Cen X-3, the 
optical depth of the VHE 7 -rays can vary from ^ 10 
to < 1 with the variation of the orbital phase from 
(/> > 0.12 (i.e. outside the pulsar eclipse) to ^ ^ 0.5 due to 
absorption on the optical/UV photons from the compan¬ 
ion. If 7 -radiation is produced close to the pulsar then a 
strong correlation of the signal strength with the orbital 
phase should be seen. 

For calculations of the orbital phase we use the 
ephemeris by Nagase et al. (1992), which are in quite 
good agreement with the ones by Ke l ley e t al. ( 1983 ) 
we used earlier (Chadwick et al. 199S , 2000 ), but being 
more contemporary are more accurate and appropriate 
for the period of our observations. Because the length of 
all ON/OFF segments accumulated during individual ob¬ 
servation nights varies and may reach several hours, i.e. 
a significant fraction of the 2.1 d orbital period, we have 
split the data of the nights with a large (> 5) overall num¬ 
ber of on-source segments into sequences of 3-4 on-source 
segments (which in some cases may overlap). In principle, 
such a grouping of data allows one to also look for varia¬ 
tions of the signal on timescales ~ 1.5 — 2 h, which could 


be of interest from a theoretical point of view (see section 
3 below). 

A trend of a uniform distribution of the 7 -ray signal 
with the orbital phase of the pulsar is apparent in Fig. 
where we show the excess events in the on-source segments 
in terms of their ratio to the number of events in the coun- 
terpa rt off -source segments. This method (see Chadwick 
et al. 1999 ) allows an assessment of the strength of a possi¬ 
ble 7 -ray signal with respect to the intensity of cosmic rays 
independent of the daily variations in the performance of 
the telescope. Note that there is some deficit in the num¬ 
ber of on-source segments corresponding to the orbital 
phase —0.5 < (p < —0.3. For all other phases, includ¬ 
ing the eclipse phase, the segments with high (say > 0 . 1 ) 
ON-OFF excess are distributed quite homogeneously. It 
is worth noting also that most of these high excess data 
points have moved further up in significance when only 
the events with the parameter brightness > 1500 digital 
counts are selected. A reasonable explanation for this ef¬ 
fect is that the 7 -ray events are indeed present in those 
segments, and that the spectrum of these 7 -rays is signif¬ 
icantly harder than the spectrum of cosmic rays. 


2.2. Pulsar spin modulation 


In our previous analyses we were looking for a modula¬ 
tion of the VHF 7 -ray signal with the X-ray pulsar in a 
very narrow band of frequencies around the fundamental 
frequency r'o = ^/Po using the entire ON-source dataset. 
Meanwhile the earlier reports on the 7 -ray signals modu¬ 
lated with the X-ray pulsar spin in Cen X-3 (Brazier et 
al. 1990a), as well as in some other X-ray binaries such as 
Cygnus X-3 (Brazier et al. 1990b) or Hercules X-1 (e.g. 
Dowthwaite et al. 1984, Lamb et al. 


had indicated 

that the pulsed 7 -ray emission was sporadic, with a du¬ 
ration typically not exceeding 1 h and occasionally signif¬ 
icantly shorter. Moreover, the frequency of 7 -pulsations 
reported in Her X-1 had on occasion been observed to be 
significantly shifted from the nominal frequency of the X- 


ray p ulsar by up to ~ 0.1 6% ( Lamb et al. 1988, Resvanis 
et al. 198g , Dingus et al. 1988 ). 


Therefore it seemed worth conducting a search for any 
possible modulation of the 7 -ray emission with the BATSE 
derived pulsar spin period in a wider bandpass around 
Pq « 4.81 s. The test initially used was the standard 
Rayleigh test for periodicity. The search was conducted 
at both the period and the half-period (in the interval 
from 4.79 s to 4.83 s and 2.395 s to 2.415 s respectively) 
due to the insensitivity of the Rayleigh test to light curves 
with a double peak separated by tt in phase; a trait occa¬ 
sionally observed in Cen X-3’s X-ray light curve (Tuohy 
1976, Nagase et al. 1992, Burderi et al. pOOO ). This feature 


suggests that both accretion hot spots of the pulsar may 
become visible during one cycle of its rotation, which is 
quite possible in the case of a relatively large inclination 
angle of the pulsar magnetic axis. It is not improbable that 
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the VHE 7 -rays might have the basic modulation near the 
half-period rather than the full spin period of the pulsar. 


The time range tested corresponds to Doppler shifted 
speeds of the putative 7 -ray source relative to the neutron 
star of up to u ^ 1200 kms“^. Note that the orbital speed 
of the Cen X-3 pulsar is 414kms“^, and the speed of the 
wind driven by the optical co mpan ion star in this binary 
is ~ 1000kms“^ (Clark et al. 1988 ). 


For a typical observation of Cen X-3 during one night 
the duration is Atobs ~ 3 hr (in a succession of ’ON’ and 
‘OFF’ segments) therefore the chosen range of trial peri¬ 
ods corresponds to ~ 18 — 19 independent Fourier fre¬ 
quency (IFF) intervals Az/ = l/A<obs around Pq (and 
about twice that for the half-period search). Given the 
wide band of trial periods studied, there is no need for an 
additional adjustment of the event arrival times to the 
pulsar position because any significant modulation fre¬ 
quency to be found could later on be easily compared 
with the Doppler shifted frequency of the X-ray pulsar at 
the known orbital position. Moreover, to adjust the event 
times to the pulsar site may well be unjustified since both 
the HE and VHE data give no indication of orbital mod¬ 
ulation, with about a quarter of all excess photons from 
the EGRET data being in the eclipse phase of the pulsar 
and the episode of highest significance in the Mark 6 data 
occuring during the X-ray eclipse. This implies the pro¬ 
duction site of the 7 -radiation to be far from the pulsar 
itself. 


We have first considered, using the standard Rayleigh 
test statistic, the data of each individual night of obser¬ 
vation of Gen X-3 after application of the full set of the 
Extensive Air Shower image parameter cuts as described 
in Chadwick et al. (1998). This analysis did not reveal any 
Rayleigh power that could be potentially significant either 
around the main period Pq ~ 4.8 s or around the second 
harmonic (P ~ 2.4s). This, however, is a result that one 
could reasonably expect. 

The Rayleigh power at some frequency v is defined 
as Pri = (c^ -I- s^)/N, where c = and s = 

sin 0 j, and 9j = 2Tri/tj is the phase of the j-ih. event. 
If the total number of events N is large and the dura¬ 
tion of the data set is the test period 1 /z^, then the 
sums c and s are practically independent (only weakly 
correlated), so for simplicity of further discussion we can 
take Pri ~ 2c^/A. For uniformly distributed events the 
mean c (and s) is equal to 0. However the root mean 
square (rms, the dispersion) of c is not zero, but is equal 
to d = = \/N/2 [|, so the expected mean Pri = 1. If 

now there are some Nq events -rays’ , G) modulated 
at frequency v, and the rest Nb = N — Nq represents 
the uniform background (cosmic rays, B) and N is the 
total number of events, then the expected mean Rayleigh 
power exceeds 1. Taking into account that the mean val- 


^ It means that c and s are expected to be in the range ±d. 


ues c ^ s ^ G/2 (for sinusoidal modulations) and that 
G ^ P ~ A, one finds: 




Ng + Nb 


0 > 


where cto = Ng/^/2.Nb- Note that, the latter variable 
represents the statistical significance of the signal in the 
data, which in practice is calculated as cr = 
which implies Aon = Nb + Ng and Aqs = Nb- 

Any cuts made will affect the number of both 7 -ray 
events, N^^ = FgNg, and background cosmic ray events, 
= FbNb, where Fg and Pr are the respective frac¬ 
tions of events surviving cuts. Since Fg can be Pr 
imaging is a very useful tool for background rejection. The 


significance of a signal, however, is cr oc P^ 


- 1 / 2 . 


therefore 


the cuts made need to be harsh to get the best signal 
to noise ratio, inevitably leading to a loss of 7 -ray events 
from the dataset. In a d.c. search this is acceptable since 
merely increasing the observation time will lead to an in¬ 
crease in the number of 7 -rays in a dataset; however, in 
a periodicity search - especially if the episodes of pulsed 
emission are short term and particularly when combined 
with a low flux of VHE 7 -rays - any loss of 7 -ray signal 
could prove fatal to a positive detection. Whilst cutting 
should improve the Rayleigh power mathematically, the 
low number of events surviving cuts means the dispersion 
of the Rayleigh power about the mean would give no reli¬ 
able estimate of the true Rayleigh power, being typically 
within the range Pri ~ Pg + ctq ± 2cjq\/Pb', where the 
sign of the last term depends on whether cb and cg are 
‘coherent’ (have the same or different signs). In order to 
get a more robust statistic the cuts need to be relaxed, to 
allow for the statistical fluctuations of Nb and Ng- 

As an example, taking the average ‘on-source’ observa¬ 
tion time to be 1.5 hours, and using even the mean com¬ 
puted flux of 1997, P(> 400 GeV) = 5 x 10“^^ cm“^s“^, 
which is a factor of 2 higher than the mean during 1997- 
1999 from Gen X-3, the mean number of 7 -rays for a 
detector with a collection area of 10® cm^ and 50% ef¬ 
ficiency (Chadwick et al. 2000| ) is ~130. Monte Carlo sim¬ 
ulations indicate that the fraction of 7 -rays passing full 
image parameter cuts is ~20%, so only 25-30 of these 
could be pulsed 7 -rays. Note that this is significantly less 
than minimum 40 — 50 which is usually needed, given the 
impact of statistical fluctuations, for a successful prac¬ 
tical application of the Rayleigh statistics, even if we 
neglect the background events. Meanwhile, if one takes 
into account also that there are on average about 400 
events remaining in each night’s data after the appli¬ 
cation of imaging cuts, the maximum expected possible 
Rayleigh power (assuming Cq + Sq — Nq) would only be 
~ 1 -I- (25)®/400 ~ 2.6. This would correspond to a 
chance probability of only ea;p(— Pri) ~ 0.07 - before the 
number of trials has been taken into account. Gompare 
this to the maximum Rayleigh power achieved by leaving 
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130 7 -ray events in a background signal of ~ 1000 events, 
found by relaxing the image parameter cuts, of 17.9 that 
would give a chance probability of ^ 10 

Thus, the hard image cut parameters which maximize 
the DC signal become destructive for the search of y-ray 
modulations with the pulsar spin period if the data of in¬ 
dividual observational nights are to be analyzed. For this 
purpose the priority should be given to keeping in the data 
as many y-ray events as possible, but trying at the same 
time to reasonably reduce the CR background events. In 
order to do so, we have applied in succession a set of signif¬ 
icantly softer image parameter cuts, removing at the first 
step from the initial data only the events with the orienta¬ 
tion parameter a > 45° (to compare with the a < 30° kept 
in the fully-cut sets). This procedure reduces by a factor of 
~ 3 the number of CR induced events, but does not signif¬ 
icantly affect the y-ray events. Analyses of these data did 
not reveal any statistically significant (i.e. corresponding 
to chance probability p below e.g. 10“^) Rayleigh peak ei¬ 
ther around the first or the second harmonic of the X-ray 
pulsar. 

The next set of data has been prepared applying a set 
of ‘soft’ image parameter cuts: by confining the position of 
the image within the camera (parameter distance)] and by 
discriminating the images using the width and the eccen¬ 
tricity parameters. This procedure can be rather efficient 
for suppressing CR background by a factor of ~ 8 — 10, 
although it may also remove some y-ray events. Lastly, 
we have also searched for pulsar periodicity signs in the 
events after discriminating the images by their brightness, 
choosing only those with > 1500 digital counts as com¬ 
pared to those with brightness > 800 digital counts. This 
procedure effectively increases the energy threshold by a 
factor ~ 2, therefore it is efficient for the timing analysis 
only if the y-ray spectrum is much harder than the dif¬ 
ferential spectrum of the local CRs with the power-law 
index acr — 2.7. This idea stems from earlier theoretical 
predictions that the spectra of episodic y-rays expected in 
the VHE domain from X-ray binaries may be anomalously 
hard due to very significant absorption on the thermal op- 
tical/UV photons produced either by the compact y-ray 
source (Aharonian & Atoyan 1991 , 1996 ) or by the op¬ 
tical companion star in the particular case of Cen X-3 
(Bednarek 2000[ ). 

The detailed analysis of all these data sets has re¬ 
vealed a strong Rayleigh power peak, corresponding to 
the chance probability for an artifact modulation of only 
p ^ 6.8 X 10“^ (see Fig. ||), in the data of observations 
on 21 Feb, 1999 at Pi = 2.3999 s. This period is blue- 
shifted from the nominal second harmonic of the X-ray 
pulsar (Po/2 = 2.408785 s, after correcting to the pulsar 
orbital phase (j) ~ 0.765 for that data) by 6.6 msec, and 
implies a motion of the y-ray source (if true) with respect 
to the neutron star with a velocity > 10^ kms“^. Interest¬ 
ingly enough, but not surprisingly, this peak is found in the 
data after both the soft image cuts and brightness > 1500 



Period (ms) 



Period (ms) 



Period (ms) 


Fig. 2. Rayleigh power probabilities in the search window 
around the half-period (the second harmonic) Po /2 of the 
X-ray pulsar in the data sets of February 21, 1999 ob¬ 
servations after application of image parameter cuts with 
different hardness. 


digital counts criteria have been applied. For comparison, 
in Fig. H we show the results for the low energy threshold 
case, brightness > 800 digital counts. The peaks disap¬ 
pear in the fully cut data, which is not surprising in light 
of the discussion above. 

A simple estimate for an overall probability ptot to find 
a Rayleigh power FRay such that po = ea;p(—PRay) < 6.8 x 
10“^ by chance can be derived taking into account that for 
a mean duration of our observations of Cen X-3 per night 
is about 3 hr, there are on average about 18.6 IFF intervals 
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Table 1. The numbers of events in the ‘ON’ and ‘OFF’ pairs of segments, the significance and the relative excess (in 
terms of ‘gamma-rays per cosmic ray’) of the signal in the data of 23 individual nights of Cen X-3 observations after 
application of the full set of the shower image parameter cuts, and for different brightnesses of the images. 



brightness > 800 digital counts 

brightness > 1500 digital counts 

brightness > 2000 digital connts 

Date 

on 

off 

sig 

gr/cr 

on 

off 

sig 

gr/cr 

on 

off 

sig 

gr/cr 

01/03/97 

446 

401 

1.546 

0.112 

265 

218 

2.139 

0.216 

199 

173 

1.348 

0.150 

03/03/97 

572 

515 

1.729 

0.111 

391 

329 

2.311 

0.188 

339 

282 

2.287 

0.202 

04/03/97 

383 

319 

2.416 

0.201 

233 

187 

2.245 

0.246 

184 

149 

1.918 

0.235 

01/06/97 

265 

256 

0.394 

0.035 

184 

171 

0.690 

0.076 

171 

156 

0.830 

0.096 

02/06/97 

455 

383 

2.487 

0.188 

304 

262 

1.765 

0.160 

260 

219 

1.873 

0.187 

04/06/97 

200 

185 

0.764 

0.081 

146 

143 

0.176 

0.021 

131 

138 

-0.427 

-0.051 

05/06/97 

323 

295 

1.126 

0.095 

218 

203 

0.731 

0.074 

182 

169 

0.694 

0.077 

07/06/97 

394 

377 

0.612 

0.045 

273 

258 

0.651 

0.058 

231 

215 

0.758 

0.074 

1997 Total 

3038 

2731 

4.042 

0.112 

2014 

1771 

3.950 

0.137 

1697 

1501 

3.466 

0.131 

27/03/98 

689 

628 

1.681 

0.097 

411 

378 

1.175 

0.087 

296 

280 

0.667 

0.057 

29/03/98 

661 

589 

2.036 

0.122 

424 

347 

2.773 

0.222 

309 

243 

2.809 

0.272 

30/03/98 

364 

373 

-0.332 

-0.024 

226 

220 

0.284 

0.027 

174 

163 

0.599 

0.067 

17/04/98 

182 

178 

0.211 

0.022 

109 

98 

0.765 

0.112 

78 

78 

0 

0 

19/04/98 

339 

336 

0.115 

0.009 

232 

216 

0.756 

0.074 

183 

167 

0.855 

0.096 

26/04/98 

59 

56 

0.280 

0.054 

45 

39 

0.655 

0.154 

39 

32 

0.831 

0.219 

27/04/98 

473 

441 

1.058 

0.073 

371 

334 

1.394 

0.111 

313 

270 

1.781 

0.159 

28/04/98 

272 

293 

-0.883 

-0.072 

206 

220 

-0.678 

-0.064 

172 

181 

-0.479 

-0.050 

29/04/98 

151 

126 

1.502 

0.198 

113 

96 

1.176 

0.177 

87 

85 

0.152 

0.024 

1998 Total 

3190 

3020 

2.157 

0.056 

2137 

1948 

2.957 

0.097 

1651 

1499 

2.708 

0.101 

13/02/99 

431 

451 

-0.673 

-0.044 

323 

337 

-0.545 

-0.042 

263 

266 

-0.130 

-0.011 

15/02/99 

78 

66 

1.00 

0.182 

42 

30 

1.414 

0.400 

30 

16 

2.064 

0.875 

16/02/99 

923 

902 

0.492 

0.023 

614 

602 

0.344 

0.020 

492 

495 

-0.095 

-0.006 

17/02/99 

922 

889 

0.775 

0.037 

643 

611 

0.904 

0.052 

536 

479 

1.789 

0.119 

20/02/99 

206 

212 

-0.293 

-0.028 

164 

170 

-0.328 

-0.035 

144 

146 

-0.117 

-0.014 

21/02/99 

355 

294 

2.39 

0.207 

248 

210 

1.778 

0.181 

210 

172 

1.944 

0.221 

1999 Total 

2915 

2814 

1.334 

0.036 

2034 

1960 

1.171 

0.0378 

1675 

1574 

1.772 

0.064 

Grand Total 

9143 

8565 

4.344 

0.067 

6185 

5679 

4.646 

0.089 

5023 

4574 

4.583 

0.098 


in the search window 4.79 — 4.83 s and twice that number 
in the second harmonic search window, 2.395 — 2.415 s. 
Because we have prepared 6 = 3x2 combinations of data 
(for different parameter cuts and brightnesses) for each of 
the 23 days, the overall number of trials is A^iff < 8 x 10^. 
This means that the overall probability of this peak arising 
by chance increases to ptot — A^iff xpo<5.4xl0“^. 

Although the chance probability is still small it is not 
a conclusive argument for VHE y-ray periodicity. In an 


attempt to gain better control of the hypothesis testing 
the analysis was repeated using a Bayes ian technique (see 
Gregory & Loredo 1992 Orford 2000 for details). This 
method applies Bayesian probability theory by comparing 
the phase distribution of a constant model for the signal 
to members of a class of models with periodic structure. 
The periodic models describe the signal plus background 
with a stepwise function, resembling a histogram, of m 
phase bins per period. In the case of a periodic model. 
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Period (ms) 

Fig. 3. Same as in Fig. |^, but for the events with param¬ 
eter brightness > 800 digital counts. 


the non-uniformity in phase is characterised by the vary¬ 
ing contents of the phase bins. Although the number of 
phase bins needed to detect any light curve and the origin 
of phase are unknowns, these ’nuisance’ parameters can 
be marginalised (integrated out); an important technical 
advantage of Bayesian inference over standard frequentist 
statistics. Since an arbitrary postulated light curve may be 
of any shape the method automatically applies Ockham’s 
razor, in that models with fewer variables are automati¬ 
cally favoured unless the evidence from the data more than 
compensates for models of greater complexity. If there are 
m phase bins the average rate A = ^ O 
fraction of the total rate in phase bin j is fj = The 
likelihood function is shown to reduce to 

( m 

n/;' 

where w is the postulated angular frequency, (j) the starting 
phase, f the set of m values of fj and Uj being the number 
of events occurring in bin j. 

The joint prior density for the parameters uj, (j), A, f is 

p{u!, 4>, A, f\Mjn) = p{uj\Mm)p{(l)\Mjn)piA\Mjn)p{i\Mjn)- 

The prior densities are 

— p{(j)\Mm) = I/Stt, this assumes any starting phase is 
equally likely. 


brightness > 1500: 

all events with a < 45° - 



395 2.4 

2.405 2.41 2.41 

Period (s) 

' 1 

soft cuts applied 

Aki 

|i , : 

395 2.4 

2.405 2.41 2.41 

Period (s) 

- 

full image cuts, a < 30° " 




2.395 2.4 2.405 2.41 

Period (s) 


2.415 


Fig. 4. Similar to Fig. g, but for the Bayesian odds anal¬ 
yses. 


— p{A\Mjn) = iMmaa:, this assumes that A does not 
change during the observation and any value of A from 
A = Q to A = Amax is possible, 

— p{uj\Mm) = where [uJhi,ujioJ\ is the least 

informative prior for the range of a;, 

— p{i\Mm) = {m — 1)!(5(1 — fj)^ where 5 denotes 
the dirac (5-function. 

The assignment of the priors beforehand is all that is 
needed before comparing the likelihoods of the models. 
The two models are equally likely a priori and so the prior 
likelihood of the non-periodic model (Mi) is p{Mi) = 1/2 
and that for a periodic model {Mm^rn = 2,mmax) is 
p{Mm\I) = l/2i^, where n = rumax - 1- 
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The final result for the odds O in favour of a periodic 
model and against a uniform model of phase when the 
phase and period are unknown (as in our case) is 


Oml — 


1 


A^!(m- 1)! r*** dio 


{N -I- m — 1)! 


U! 


n2T: 

0 


that would first of all secure a larger number of possible 
7 -ray events in the remaining data, could possibly have a 
better chance to enhance the initial signal. 

3. Possible models 


where 4>) is the number of ways that the set of nj 

counts can be made by distributing N counts in m bins: 




Nl 


nr=i^.! 


and rij, the number of events placed in the jth phase bin 
depends on uj, (j) and m. The overall likelihood of period¬ 
icity is then found from 

rrimax 

Oper — ^ ^ OjYi\ . 

m— 2 

The attractiveness of the Bayesian technique is that 
the overall probability of the hypothesis of a time- 
modulated signal does not scale linearly with the number 
of trials as in the Rayleigh test, but with the ln(period 
range searched). The periodicity search is therefore only 
penalised by the fact that the period is uncertain within 
a range and not by the number of times that region is 
searched. In Fig. ^ we show the Bayesian odds results for 
the same data sets as used in Fig. for the Rayleigh statis¬ 
tics. Formally, the peak in the median panel in Fig. ^ 
means that the hypothesis that the data set does con¬ 
tain a signal modulated at the period Pi = 2.3998 s is 
.^odds = 1.2 X 10^ times more probable than the hypothe¬ 
sis of a uniform time distribution. Note that the positions 
and the qualitative behaviour of the 2 highest peaks in 
Figures || and ^ are similar. 

Finally, it is worth noting the earlier reports in 
the 1980’s about the detection of anomalously pulsating 
episodic VHE y-ray emission from Her X-1 by three in¬ 
dependent groups (Whipple, Lamb et al. 1988; Haleakala, 
Resvanis et al. 1988; and Los Alamos ,Dingus et al. 1988), 
offset from the main period by ~ 0.16%. It is also worth 
noting that Her X-1 had a short episode (At < Ihr) of 
y-ray emission reported as being detected simultaneously 
and independently by the University of Durham (Chad¬ 
wick et al. 1987) and the Whipple Observatory (Gorham 
et al. 1986 ) groups, with both pulsation frequencies and 
estimated fluxes being consistent in the initial uncut data 
of both groups. After application of hard AZWIDTH pa¬ 
rameter cuts which had earlier produced the best results 
for the Crab Nebula signal (Weekes et al. 1989), the pulsa¬ 
tions in all Her X-1 data of the Whipple group disappeared 


(Reynolds et al. 1991). This, however, is just the outcome 
that one would expect, as we have discussed above in this 
Section. Taking into account the rather steep spectrum of 
the Crab Nebula compared to the one expected from Her 
X-1 (Aharonian & Atoyan 1996 ), we may here hypoth¬ 
esise that an application of much softer parameter cuts. 


3.1. General considerations 

The results of observations of Cen X-3 in both HE and 
VHE y-ray domains allow us to significantly constrain dif¬ 
ferent models of y-ray production in this X-ray binary. 
The first informative result is the absence of any signs of 
orbital modulation of the y-ray signal in either of these en¬ 
ergy domains. Because of the presence of a very luminous 
0-type star close to the pulsar the VHE radiation pro¬ 
duced at the pulsar orbit would be significantly absorbed 
on the optical/UV photons from the companion prior to 
escape from the system. Therefore the y-ray source can¬ 
not be placed close to the neutron star. This conclusion is 
also in agreement with the fact of detection of y-rays by 
EGRET (Vestrand et al. 1997), and probably also by the 
Mark -6 telescope, in the phase of complete eclipse of the 
X-ray pulsar. 

Thus, at least for this particular binary any model 
assuming production of y-rays directly by the pulsar, or 
models invoking the accretion disc around the neutron star 
as a s ite for y-ray pr oduct ion (e.g. Chanmugam & B reche r 
1985, Katz & Smith 1988, Raubenheimer and Smit 1997), 
can be excluded from the list of reasonable options. At 
the same time, the y-ray source must be genetically con¬ 
nected to the neutron star. This requirement follows not 
from the effect of y-ray modulations with the pulsar spin, 
which would still need additional confirmation by future 
observations, but from the general considerations of the 
source energetics. 

Indeed, the radiation flux detected by EGRET during 
2 weeks in 1994 corresponds to an average source lumi¬ 
nosity of L.y(100MeV < E < lOGeV) ~ 5 x 10^®ergs“^ 
(Vestrand et al. 1997), and the mean flux detected in 1997- 
1999 by the Mark -6 telescope corresponds to a luminosity 
L.y{> 400GeV) ~ 10^®ergs“^. For the parent relativistic 
particles, either electrons or protons, these luminosities in¬ 
evitably imply an acceleration power Pace 10 ^^ ergs“^ or 
higher. The optical/UV luminosity of the Krzeminski star 
corresponding to the b olometric magnitude Afboi — —9 
(Hutchings et al. 1979|) is really high, Lboi — 10^® er g/€l 
It would suffice if a small fraction of this luminosity could 
somehow be converted to the acceleration of relativistic 
particles. In principle, this might be due to a radiatively 
driven supersonic stellar wind and subsequent production 
of shocks in such winds. Meanwhile, for the characteristic 
speed of the stellar wind v ~ 10 ^ kms ~^ and the mass-loss 
rate M ^ lO“®M 0 yr“^ (Glark et al. 1988) the entire ki- 


^ Note that for a massive star with M ~ 20 Mq this is al¬ 
ready comparable with, but still significantly below the Ed¬ 


dington luminosity Lsdd = 1.3 x 1O'^“(M/M0) ergs 


„-i 
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netic energy of the wind produced by the Krzeminski star 
makes only Pwind = Mv‘^/2 ~ 3 x 10^^ergs“^. 


Any reasonable possibility for production of a kinetic 
power exceeding 10^^ ergs“^ in Cen X-3/V779 Cen binary 
to be is therefore to be attributed to the neutron star. 
For the 4.8 s spin period of Cen X-3 the classical mecha¬ 
nism of the magnetic dipole radiation falls several orders of 
magnitude below the values required. Therefore the only 
remaining principal option for the prime energy source for 
the acceleration of relativistic particles in this system re¬ 
mains the kinetic energy of the inner accretion disc formed 
around the neutron star. Most of the kinetic energy of the 
disc ends up on the neutron star surface in the form of 
the thermal plasma energy responsible for the observed 
X-rays with the luminosity reaching (in the high state) 
Ax > 10^®ergs“^. It is possible, however, that a signifi¬ 
cant fraction of the accretion disc energy can be ejected 
from the system producing powerful outflows in the form 
of two-sided jets observed at least in some class of close 
X-ray binaries currently identified as ‘microquasars’ (see 
Mirabel & Rodriguez 1999 for a review), which may con¬ 
tain either a stellar-mass black hole or a neutron star. 
A spectacular example of a microquasar containing an ac¬ 
creting neutron star that produces powerful subrelativistic 
(ujet ~ 0.26c) jets is SS 433, with the jet kinetic energy es¬ 
timated in the range from ~1^® ergs“^ to 4 x 10 ^° ergs“^ 
(see Margon & Anderson 1989 , Panferov 1999 ). In this 
regard it is also worth noting the recent report by Jerni- 
gan et al. ( 2000 ) of the detection of high-frequency QPOs 
(‘quasi-periodic oscillations’) in the X-ray emission of Cen 
X-3, which is a characteristic feature of microquasars. 


The principal scenarios for y-ray production in Cen X- 
3 are thus connected with powerful jets of plasma driven 
by the inner accretion disc. In one such scenario these jets 
would create strong shocks while propagating through the 
rather dense wind driven by the 0-star. Subsequent ac¬ 
celeration of particles on these shocks could result in the 
creation of a y-ray source around the region of jet prop¬ 
agation (or its damping) on large spatial scales Is compa¬ 
rable with, or more probably significantly, exceeding the 
size of the binary ~ 10^^ cm. Obviously, such a spatially 
extended model will not be able to explain episodes of 
pulsed y-ray emission that may have been detected from 
the system. Evidence for pulsed emission from Cen X-3 
is far from conclusive, but if it is there it probably con¬ 
stitutes only a fraction of the overall y-radiation of Cen 
X-3. An interpretation of the ‘anomalously’ pulsed y-ray 
episodes may be possible in the framework of compact- 
source model scenarios that would suppose direct interac¬ 
tion of the jet either with material ejected from the inner 
parts of the accretion disc, or with a dense external target 
that might accidentally fall under the jet. Below we will 
discuss these model possibilities. 


3.2. Spatially extended source model 

In the framework of this model that assumes acceleration 
of particles on the shocks driven into the supersonic (and 
therefore favourable for generation of shocks) wind of the 
0-star by jets appearing from the inner accretion disc, 
the source size 4 ^ 10^^ cm, so that the y-ray source is 
surrounding the X-ray pulsar. This circumstance allows us 
to exclude with a significant confidence the hadronic (tt^- 
decay) origin of y-rays for the spatially extended model. 

Indeed, the energy loss time of protons due to inelastic 
collisions with ambient gas with a density rigas (in terms 
of nucleons, or ‘H-atoms’) is 

( Ti \ — ^ 

IQio cm“3/ ^ 

where Kp « 0.45 is the inelasticity coefficient and CTpp ^ 
30 mb is the cross-section of nuclear interactions of protons 
at GeV energies which slowly increases with the particle 
energy. Taking then into account also that in pp-collisions 
only ~ 1/3 of the initial energy of relativistic particles 
is eventually transformed to y-rays (the rest goes to the 
secondary electrons and neutrinos), the total energy Wp of 
relativistic protons which is needed for y-ray production 
with luminosity A-y is estimated as 

Wp ~ 3tpp X Lj (2) 


For a y-ray source extending up to characteristic 
distances Rg measured from the centre of the 0-star, 
the source volume can be generally written as 14 = 
/v X 47ri?g/3 where fy < 1 is the volume filling factor. 
The energy density of relativistic protons Wp = Wp/Vg 
should then be compared with the energy density Wg^ = 
GMrr^rig^g thermal gas in the gravitational field of 

the massive companion star with M ~ 20Mq which de¬ 
fines the gravitational confinement potential of the binary. 
Obviously, in any reasonable scenario the energy density 
of relativistic charged particles that can accumulate {be 
frozen) in the ambient plasma should not exceed the gravi¬ 
tational field energy density, otherwise the relativistic par¬ 
ticle gas cannot be confined to a scale Rg but would act to 
inflate the source, also strongly affecting the gas dynamics 
in the binary system. Then from the condition Wp < Wgr 
we can estimate the maximum y-ray luminosity AJ/o^ that 
can be produced by relativistic protons: 


47r o 

—fyGMmpKpappN^c 

y 


(3) 


~ 2.5 X 10^^ fy 


M 


20 Mp 


iVn 


10^^ cm“^ 


ergs 


-1 


Here A^h = ?T-gas x Rg is the gas column density of the y- 
ray source which we have normalized to 10^^ cm“^. This is 
a reasonable maximum value for A^h that corresponds to 
the gas density rigas < 10^^ cm~^ of the wind at distances 
Rs > 10^^ cm (Clark et al. 1988 ). Eq.(3) shows that the 
maximum luminosity to be expected in y-rays of hadronic 
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origin is orders of magnitude below the 7 -ray luminosities 
detected from Cen X-3. 

It should be noted that Eq.(3) is actually a rather uni¬ 
versal upper limit for the 7r°-decay 7 -ray luminosity that 
can be produced {quasi-stationarily) by a source gravi¬ 
tationally confined to a central mass M. In particular, 
this relation is valid for 7r°-decay 7 -rays that can be pro¬ 
duced in the accretion discs of neutron stars or black holes 
{including supermassive BH in AGNs). For Cen X-3 we 
should substitute M ~ 1.2 Mq and take into ac¬ 

count that for any reasonable disc geometry /y cannot 
exceed ~ 0.1. Then even assuming A^h — 8 x 10^^cm“^ 
which is the characteristic maximum column density accu¬ 
mulated across th e accr etion disc of Cen X-3, as deduced 
by (Nagase et al. 1992|) from the X-ray line variations in 
the pre-eclipse dips, we conclude that any hadronic model 
invoking the accretion disc as a site for 7 -ray produc¬ 
tion would be far too inefficient for an explanation of the 
observed fluxes, or otherwise the relativistic proton gas 
would blow up the accretion disc itself. 

Considering now models that assume a leptonic ori¬ 
gin of 7 -rays, it becomes crucial for the existence of a 
quasi-stationary large-scale source that the energy loss 
time of the electrons in Cen X-3 is much shorter than 
tpp. At distances R > 1.5i?o from the 0-star the en¬ 
ergy density of the thermal optical/UV radiation field is 
w\jY = Aboi/dTTtherefore the cooling time of an elec¬ 
tron with Lorenz-factor 7 due to inverse Compton (IC) 
radiation can be estimated as 


tic - 1-1 X 10 '‘(i?/ 10 ^^cm) 27 “i s 


(4) 


This estimate is valid in the Thompson limit for IC scat¬ 
tering when the parameter b = 47(eo/TOeC^) < 1. For 
the 0-star in Cen X-3 with T ~ 35 000 K the mean en¬ 
ergy of the thermal photons is cq ~ 3 kT ~ 9 eV, there¬ 
fore the transition from the Thompson limit to the rel¬ 
ativistic Klein-Nishina limit of IC scattering occurs at 
7 ~ 7kn = 1.4x 10'^. Thus, for electrons with 7 ^ 10^ —10“^ 
responsible for production of high-energy 7 -rays the IC 
cooling time is orders of magnitude smaller than the cool¬ 
ing time of relativistic protons given by Fq.(l). For elec¬ 
trons with 7 > 7 kn the IC loss time does not drop further 
as Fq.(4) would predict, but rather increases, so that for 
VHF electrons the IC cooling time at distances R is about 
(10^ — 10^) X cm)^. Thus, for these electrons as 

well the case is still tic "C tpp, which leads to a total 
energy We needed in relativistic electrons which is much 
smaller than Wp. 

Thus, a spatially extended 7 -ray source assumes dif¬ 
fusive shock acceleration (e.g. see Drury 1983, Blandford 
& Eichler 1987) of the electrons on spatial scales Rs com¬ 
parable or exceeding the characteristic length scale of the 
binary. In order to estimate a plausible size Rs, let us es¬ 
timate the characteristic time needed for acceleration up 

® Note that with increasing R the gas density rigas drops 
approximately as oc R~^, so tpp oc R^ similar to tic. 


to energies E ^ lOTeV. Using the results of Lagage & 
Cesarski 1983| , in the Bohm diffusion limit this time can 
be expressed as 


tacr, - 7 X 10 " 


E 


lOTeV 



1 V ) 

VigJ 

' V 3000 kms " J 


-2 


(5) 


where Ush is a mean shock speed. Since the characteristic 
speed of the stellar wind in Cen X-3 is Uw ~ 10^kms“", 
the normalization of Ush to ~ 3000 kms“" used above 
seems reasonable. Now taking into account that because 
of advection in the plasma wind all relativistic particles 
leave the radial scales R at least on the convective escape 
timescale tesc — R/v^, the characteristic size of the region 
required in this scenario for production of multi-TeV elec¬ 
trons should satisfy the condition tacc(lOTeV) < tesc- As 
it follows from our calculations below, the magnetic fields 
of order U ^ 1 G in the stellar wind of Cen X-3 would 
not contradict the fact of non-detection of radio fluxes 


from Cen X-3 (see Vestrand et al. 1997). Therefore, from 
the condition tgsc > Ccc, the size of the TeV 7 -ray source 
is estimated as i?s > 5 x 10"^cm. On the other hand, the 
source size should not significantly exceed 10 "" cm because 
of the fast decline of the density of the UV photon field 
from the 0 -star, which essentially reduces the efficiency of 
the IC 7 -ray production (as tic > fesc)- 

Note that Rs is much larger than the radius Rq = 
8.6 X 10"" cm of the Krzeminski star or the distance 
a = 1.3 X 10"^ cm between the stars. This implies a source 
volume filling factor /y ~ 1. More importantly, the large 
Rs suggests that the 7 -ray source is quasi-stationary on 
time scales > 1 day. 

Provided that the volume occupied by the shock fronts 
(i.e. where the particle acceleration occurs) is significantly 
smaller than the 7 -ray source entire volume, one can 
describe the energy distribution of relativistic electrons 
N{E) by a kinetic equation (see e.g. Ginzburg & Sy- 
rovatskii 1964): 


at dE^ > tesc 


Q, 


( 6 ) 


where Q = Q{E,t) is the rate of the electron accelera¬ 
tion/injection, and P = —dE/dt corresponds to overall 
electron energy losses averaged over the volume of the 
source. The injection spectrum of the electrons is approx¬ 
imated as 


Q{E) cx exp {-E/Eo) , (7) 

with Og ~ 2 typical for the diffusive acceleration on the 
strong shock fronts, and Eq > 10 TeV for the exponential 
cutoff energy. 

In Figure we show the spectra of non-thermal radi¬ 
ation from radio wavelengths up to the HF 7 -ray band 
calculated for such a quasi-stationary spatially extended 
source with i?s = 6 x 10 "^ s, assuming injection of relativis¬ 
tic electrons with an overall power Pace = 2 x 10^"" erg. In 
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Fig. 5. The spectra of synchrotron (solid lines), 
bremsstrahlung (dot-dashed lines), and inverse Comp¬ 
ton (dashed lines) radiations expected when calculated 
in the framework of a quasi-stationary model with Rs = 
6 X 10^^ s, assuming a magnetic field B = 0.7 G, and in¬ 
jection spectrum of electrons with ag = 2.0, Eq = 15 TeV, 
and Pace = 2 X 10^^ erg. Radiation fluxes produced outside 
the source (zone 2, see text) are shown by thin lines. 


order to assess a contribution to the overall fluxes due 
to electrons taken away to large distances with the stel¬ 
lar wind, the calculations are done in a two-zone model 
approach (see Atoyan et al. 2000 for details), which as¬ 
sumes that there is no acceleration of fresh electrons at 
R > Rs (zone 2), but there is rather an injection of elec¬ 
trons escaping from the main source (zone 1) at the rate 
Q' = N{E)/tesc- In Figure ^ the fluxes of synchrotron, 
bremsstrahlung and IC radiation produced in zone 2 are 
shown by thin lines. Obviously these fluxes are negligi¬ 
ble compared with the respective ones produced in the 
region R < Rg where the effective electron acceleration 
is supposed. It tells us that for most of the electrons the 
characteristic cooling time inside the source is indeed sig¬ 
nificantly shorter than their escape time. This is apparent 
also in Figure ^ where the energy distribution of relativis¬ 
tic electrons in zone 2 is shown by a dot-dashed line. Note 
an unusual profoundly concave shape of the electron dis¬ 
tribution which is due to the Klein-Nishina effect in the 
IC energy losses of VHE electrons in the thermal radiation 
held of the 0-star. 

The full dot in Figure ^ corresponds to the upper flux 
limit Sn < 70m Jy at 5 GHz in the direction of Cen X-3 
(Vestrand et al. 1997 ). The heavy solid curve in this fig¬ 
ure shows that due to synchrotron self-absorption at fre¬ 
quencies below 10-15 GHz the values of the magnetic field 
up to ^ 1 G at distances < 10^^ cm can still be accepted. 
However an assumption of a stronger magnetic held would 
result in significant synchrotron fluxes at sub-mm wave¬ 
lengths which would then already be observed. 

In Figure ^ we also show a typical level of X-ray fluxes 
detected by BATSE from Cen X-3 in the active (‘high’) 
state. These fluxes are predominantly of thermal origin, 
and much higher than the non-thermal fluxes. Note how¬ 



Fig. 6. The energy spectra of relativistic electrons formed 
in zones 1 (dashed) and zone 2 (dot-dashed), and the total 
spectrum (solid line) for the extended source model with 
parameters in Fig. || 


ever that the latter could, in principle, show up in the 
hard X-ray/soft y-ray band, E ~ (0.1 — 10) MeV, where 
the thermal radiation of the neutron star drops. In the 
region of HE y-rays the IC radiation flux a ~ 2 is some¬ 
what steeper than the October 1994 flux of Cen X-3 with 
a hard mean differential power-law index aobs = 1-81 de¬ 
tected by EGRET (hatched region in Figure ||). However, 
if we take into account the reported uncertainty in the 
power-law index ±Aa = 0.37, the model appears in good 
agreement with the EGRET data (see Figure |^, where 
the upper boundary of the hatched zone at 0.1-10 GeV 
energies shows the median flux, while the lower one cor¬ 
responds to ttobs = 2.18). 

In Fig. ^ IC y-ray fluxes extending beyond 10 TeV 
are presented. The solid curve corresponds to the unab¬ 
sorbed flux produced in the source. The principal tar¬ 
get photons for the IC scattering of the electrons in this 
model is the thermal UV radiation of the 0-star, the con¬ 
tributions due to X-ray (3-dot-dashed) and synchrotron 
(dot-dashed) photons being negligible at all E. The heavy 
dashed line shows the fluxes of y-rays escaping the source, 
after taking into account their absorption in the UV pho¬ 
ton held. Calculations assume that y-rays are produced 
uniformly at the scales R < Rg. 

In order to demonstrate a possible impact of pho¬ 
toabsorption, in Figure we show by full dots the flux 
that would escape the binary if the same unabsorbed y- 
radiation were to be produced close to the neutron star 
in the orbital phase (p ^ ±0.25. In that case the optical 
depth Tyj{E) would reach a maximum ~ 7 at energies 
E ^ 200 GeV, and yy-absorption would remain very sig¬ 
nificant in the entire region from E > 30GeV. Because 
depends strongly on p, essentially becoming less than 1 at 
a phase of ~ 0.5, a strong orbital modulation of VHE 
y-ray fluxes produced close to the neutron star would be 
expected. In our calculations we did not take into account 
a possible contribution connected with the development 
of pair-photon cascades in the atmosphere of Cen X-3, 
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Fig. 7. The spectra of IC 7 -rays calculated for the spa¬ 
tially extended source with parameters in Fig. |^. Contribu¬ 
tions to the total unabsorbed IC radiation (solid line) due 
to upscattering of the UV radiation of the companion O- 
star (thin dashed line; distinguished only at E < 10® eV), 
X-ray pulsar photons (3-dot-dashed line) and synchrotron 
photons (dot-dashed line) are shown. The heavy dashed 
line shows the spectrum escaping the source. The heavy 
dotted line shows the fluxes that would be expected if the 
same unabsorbed radiation were produced at the pulsar 
orbit; the stars show the unabsorbed radiation of hadronic 
origin (see text). The hatched region at TeV energies cor¬ 
responds to the average flux detected by the Mark -6 tele¬ 
scope during 1997 (which is higher than the mean for 1997- 
1999, but closer to 1994) calculated for differential power- 
law indices between = 2 and 3. The range of differential 
fluxes, from median to low, detected by EGRET during 
October 1994, as well as the upper flux limit for earlier 
observation period, are also shown. 


which in practice remains small and cannot in any way 
significantly compensate the impact o f abs orption on the 
spectra of VHE 7 -rays (see Bednarek 2000 ). 


The total amount of energy accumulated in relativis¬ 
tic electrons shown in Figure || is We = 2.4 x 10'^° erg. 
For comparison of the efficiencies of hadronic and lep- 
tonic models, in Figure ^ we show by stars the fluxes of 
7r°-decay 7 -rays produced by relativistic protons with the 
same amount of total energy, Wp = 2.4 x 10'^° erg. 


It should nevertheless be said that the energy We is 
still significantly larger than the total potential energy 
Wq of the plasma in the gravitational field of the binary 
at scales < 10^® cm. The latter can be expressed directly 
through the mass loss rate M = WmpUgasR^Vvj of the O- 
star if we approximate the stellar wind speed as a constant, 
Vv,{R) — 10 ®kms“^: 




(8) 


where M-q = M/10 ^ MQ-yi However, We is compa¬ 
rable with the total kinetic energy of the wind 


Wkin ^ 3.2 X 10^° M_6 


R 


10 ®kms“i 10 ^® cm 


erg. 


(9) 


It then implies that this model can still be magnetohy- 
drodynamically self-consistent because the pressure of the 
relativistic electron gas can be ‘balanced’ by the kinetic 
(ram) pressure. Obviously, the electrons in these circum¬ 
stances should contribute very signihcantly to the gas dy¬ 
namics of the stellar wind. 

It is worth comparing here requirements that would 
follow from similar arguments in the case of a spatially 
extended hadronic model. In order to explain the obser¬ 
vations one would require an overall energy in relativistic 
protons almost 3 orders of magnitude larger than the flux 
of gamma-rays of hadronic origin (shown by stars in Fig. 
^); i.e. Wp would exceed 10'^®erg. Such a huge amount of 
energy at a scale of Rg =< 10^® cm, as assumed for Fig. ^ 
is very problematic to expect from a binary star generally, 
and is far beyond the gravitational confinement abilities 
of the star, as discussed above. Moreover, Eq.(9) predicts 
that one would then expect such a relativistic proton gas 
would drive a wind in Cen. X-3 moving at subrelativistic 
speeds (z^w)- Assuming a larger source size would not help, 
as might be hinted by the same Eq.(9), on the contrary it 
will make the situation even worse because of a fast drop 
off of the mean gas density at larger distances resulting in 
a fast increase of the overall energy of relativistic protons 
needed. Only a rather compact source that could provide 
gas densities much higher than are found in the 0 -star 
wind, so resulting in a reasonable Wp, giving a reasonable 
chance to provide the observed fluxes (see below). Note, 
however, that this does not mean that in the framework 
of a large-scale model protons are not allowed to be accel¬ 
erated with an efficiency comparable with the electrons. 

The relativistic electron gas pressure in the source can 
be reduced further if we take into account that most of the 
energy in Figure ^ is due to particles with E < 100 MeV 
which do not contribute to production of HE 7 -rays. Then 
we can reduce We assuming a very hard spectrum of ac¬ 
celerated electrons with ae < 2. Note that taking into ac¬ 
count non-linear modification effects at the strong shock 
fronts, power-law spectral indices of the shock acceler¬ 
ated particles even as hard as ae = 1.5 could be expected 


(Maikov 1997). 

In Fig. 1^ we show the 7 -ray spectra calculated for 
ae — 1.75, and assuming an injection rate of electrons 
with Pace = 6 X 10®®ergs“^. In this case the total en¬ 
ergy of the electrons in the region R < Rg = 10 ^® cm is 
only We = 4.1 x 10®® erg. Note that in the HE 7 -ray do¬ 
main the shape of the model spectrum now agrees well 
with the median spectral index of the fluxes detected by 
EGRET in October 1994. Meanwhile, the absolute fluxes 
of HE 7 -radiation are now in agreement with the upper 
flux limits reported for the 9 month earlier observing ses- 
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Fig. 8. Gamma-ray fluxes (heavy dot-dashed line) ex¬ 
pected in the framework of a quasi-stationary source with 
Rs = 10^^ cm, i? = 1 G in the case of a hard power-law 
spectrum of accelerated electrons with Oe = 1.75, and 
total electron acceleration power Pace = 6 x 10^®ergs“^. 
Gontributions due to IG (dashed line), synchrotron (thin 
dot-dashed line) and bremsstrahlung (3-dot-dashed line) 
radiation mechanisms are shown. 


sion of Gen X-3 by EGRET (Vestrand et al. 1997). For 
calculations in Fig. || we have chosen rather a small accel¬ 
eration rate of electrons. Pace = 6 x 10^®ergs“^, in order 
to demonstrate that in principle in the framework of this 
quasi-stationary model with Oe < 2 it is possible to ex¬ 
pect episodes (on time scales > 1 d) of y-radiation that 
could have been detected in the VHE y-rays, but would 
be missed by EGRET. 

In conclusion, a spatially extended quasi-stationary 
source model seems capable of providing a phenomenolog¬ 
ically self-consistent explanation for most of the currently 
existing data of y-ray observations of Gen X-3. The single 
observational feature which cannot be addressed by this 
model is an interpretation of possible modulations of the 
y-ray emission with the pulsar spin period Pq. If genuine, 
such modulations would either require an alternative or 
an additional much more compact y-ray source. 


3.3. Compact source models 

An obvious requirement on any model that would be able 
to explain a y-ray emission modulated at the pulsar spin 
period Po is that the effective source radius Rs has to 
be noticeably smaller than cPo/2. For Gen X-3 therefore 
Rs should be < 5 X 10^° cm. Another requirement is that 
the energy loss time tioss of relativistic particles produc¬ 
ing pulsed y-rays should also be smaller than Pq. Then 
for floss ^ 1 s the energy density of relativistic particles 
can be estimated as Wrei > lO^erg/cm^ using the relation 
kErei > L.ytioss for the total particle energy in the source. 
Such high energy densities can be reached in the binary 
only in the inner accretion disc around the neutron star, 
which however cannot be a site for the y-ray emission 


because of the absence of any indications for the orbital 
modulation of the fluxes in both the HE and VHE do¬ 
mains, as discussed in Section 2 above. This implies that 
the source of pulsed y-radiation in Gen X-3 cannot be con¬ 
fined, so that generally one would expect only rather short 
episodes of pulsed y-ray emission. 

As we have discussed above, a source of the observed 
energetic gamma-radiation in this X-ray binary should 
be connected with jets generated by the accretion disc 
around the neutron star. There are then two principal 
types of model, leptonic and hadronic, for the production 
of episodic pulsed radiation, both of which would assume 
a compact source (‘clouds’, ‘blobs’) propagating in the jet. 
Pulsations could be produced if a powerful relativistic en¬ 
ergy outflow, in the form of either electromagnetic Poynt- 
ing flux or a beam of relativistic particles, modulated with 
the pulsar spin period propagates in the jet region and ac¬ 
celerates/injects relativistic particles in the source. 

A hadronic ‘beam-target’ model for production of 
pulsed VHE y-radiation in X-ray binar ies h a s bee n sug¬ 
gested earlier by Aharonian & Atoyan ( |l991 , 1996 ). This 
model assumes that a powerful beam of relativistic pro¬ 
tons accelerated in the vicinity of the pulsar hits a dense 
plasma cloud that may appear in the jet propagation re¬ 
gion. As it follows from Eq.(l), the gas density in the cloud 
should be very high, rigas > 10^®cm“^, in order to pro¬ 
vide fast energy losses for the protons, tpp < Pq = 4.8 s. 
For a cloud with a radius i?ci ~ 10^° cm this implies a 
cloud mass Md > 10^^ g. One may phenomenologically 
suppose that either such dense plasma blobs are ejected 
from the massive optical companion star, and then they 
may cross a conical region of jet propagation. We may 
also speculate that compact dense plasma blobs can be 
ejected from the accretion disc or be created in the jet re¬ 
gion, and then the ‘beam-target’ interaction takes place in 
the episodes when relativistic protons are sporadically ac¬ 
celerated in the central engine (e.g. by the dynamo action, 
magnetic field line reconnection, etc) forming relativistic 
proton beams streaming in the jets, this regard that very 
dense and compact (~ 10® cm) gas clouds responsible for 
the observed optical emission in the powerful relativistic 
jets of the prominent binary SS 433, are supposed to be 
produced/condensed due to thermal instabilities in those 
jets (see Brinkmann et al. 198S , Panferov 1999|) . 

A relativistic proton beam interacting with the gas in 
the cloud would produce Tr'^-decay y-rays. Non-thermal 
radiation in this model will be contributed also by sec¬ 
ondary electrons (e^). At the same time some fraction 
(a few per cent, see Aharonian & Atoyan 1991) of the 
injected proton energy would inevitably go to Goulomb 
(‘ionization’) losses heating the cloud up to temperatures 
Tci ~ (5 — 10) X 10^ K. This results in a very high density 
UV radiation field in the cloud, so that its opacity with 
respect to VHE y-rays can be very high. The spectra of 
escaping y-rays are then essentially defined by the pair- 
photon cascade developing in the cloud. At the same time. 
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because of the high pressure created, the source would ex¬ 
pand with a speed presumably of order of the sound speed, 
^ ^kTci/irip, if not confined by some external pressure. 
Therefore the radiation spectra predicted in this expand¬ 
ing source should be rapidly evolving. 

In Fig. 1^ we show the non-thermal radiation spectra 
expected in this model at 4 different times t, from 5 min to 
1 day (= 86400 s), after a cloud with a mass Md = 10^^ g 
and an initial size i?ci = 5 x 10® cm is hit by a powerful 
beam of relativistic protons at the distance D^i = 8 x 
10^^ cm from the neutron star. The method of calculations 
used and other details of the model are found in Aharonian 
& Atoyan (1996). The spatial angular density of the beam 
power needed in this model is rather high, dPbeam/dSl = 
3 X 10"*^^ ergs“^sr“^ in Fig. This implies that in order 
to have reasonable total energetics, the relativistic proton 
beam should be collimated to within a few degrees. Thus, 
for the collimation angle 5° the total beam power would 
be Pbeam = 1-8 X 10®®ergs“^. We remind the reader in 
this regard that the collimation angle of SS 433 jets is 

< 1.4°, and the jet kinetic energy is estimated up to 
several times 10"'^® ergs“^ (e.g. Margon & Anderson 1989|) . 

Under the pressure of such a powerful beam the ex¬ 
panding clouds would be expelled from the binary reach¬ 
ing speeds of order lO'^ kms“^. This is seen in Fig. where 
we show the increasing distance Dc\ of the cloud from the 
neutron star. The characteristic proton cooling times and 
the cloud’s radius are also shown. For the cloud mass sup¬ 
posed in this case one could expect a pulsed emission only 
during the Hrst < 1 /i, since at t = 3 h the gas density in 
the cloud drops below 10^® cm“® so that tpp > Pq- Note, 
however, that calculations in Fig. assume an external 
pressure free expansion of the cloud with the sound speed. 

In fact the expansion of a fast cloud may be signif¬ 
icantly confined due to the ram pressure of the external 
medium Pram = »7ipngasU^i. For a gas density in the wind of 
the 0 -star declining as the maximum distance scale 
(.Rconf) in the binary up to which the ram pressure will 
be able to conhne the cloud’s expansion can be estimated 
from the condition Pram > Prei, where Prei = u;rei/3 is 
the pressure of relativistic particles. For iCrei > lO^ergs”^ 
estimated above this results in 


Pconf < 5 X 10^2 M. 


1/2 / 


■Ccl 


-6 


V 10 ^ kms" 


cm 


( 10 ) 


where we have assumed Uw = 10 ®kms“^ for the wind 
speed. Thus, for a mass loss rate M < 3 x lO“®M 0 yr“^ 
the source confinement time 


Uonf = < 5 X lO^Ml/' 


^'cl 


( 11 ) 


does not exceed a couple of hours, independently of the 
source speed. However, a high speed of source is impor¬ 
tant in order to provide Pconf ^ 10 ^^ cm because other¬ 
wise the flux of pulsed VHE y-rays would be effectively 
suppressed. Eq.(lO) then suggests that Ud ^ lO^kms”^ or 
larger. Note that both Eqs.(lO) and (11) limitations relate 


not only to the hadronic model but in particular to the 
leptonic model for interpretation of pulsed y-ray emission 
considered below in this section. 

An interesting feature in Fig. ^ is a profound absorp¬ 
tion of y-radiation in the thermal UV radiation field of the 
cloud itself, and not only from that of the 0-star. It leads 
to an unusually hard spectrum of VHE y-rays predicted by 
this model at all initial stages t < (1 — 2 ) h when the pulsed 
signal is produced. In this regard, it is worth remembering 
that the most significant indications for a pulsed signal in 
our 21 February 1999 data are found in the air shower 
events corresponding to enhanced brightness. 


Considering now a compact source with a leptonic ori¬ 
gin for the pulsed y-radiation, note that in this case we 
cannot assume that electrons are accelerated close to the 
pulsar and then supplied in a relativistic beam into the 
cloud because the life time of HE and VHE electrons in 
the radiation held of the companion star (Eq.4), is much 
shorter than their travel time to distances 10 ^^ cm 
where the source of VHE y-rays should be located, as dis¬ 
cussed above. 

The principal leptonic model for Cen X-3 then basi¬ 
cally corresponds to the model developed earlier for the 
microquasars, i.e. galactic sources with relativistic jets like 
GRS 1915-1-105. For a detailed de script ion of the model 
we refer to Atoyan & Aharonian ( 1999| ). Generally, this 
model assumes that the inner accretion disc of the com¬ 
pact object in the binary, in our case the neutron star (or 
a galactic BH candidate for GRS 1915-1-105) sporadically 
ejects a pair of clouds in two opposite directions to each 
other. In the case of stellar-mass black holes the speed of 
the ejecta may reach typically > 0.9c, whereas in the case 
of neutron star discs significantly smaller jet speeds are 
assumed (see Mirabel & Rodriguez 199S). These clouds 
should then be energized from the central engine by the 
relativistic wind (e.g. Poynting flux) propagating in the 
jet region and modulated with the pulsar spin period. 
The latter condition would not be needed if the pulsed 
y-rays should not be explained. Note that certain indica¬ 
tions that the pairs of ejecta in GRS 1915-1-105 are contin¬ 
uously fed with relativistic energy from the central source 
are found in the data of radio observations of this micro¬ 
quasar (see Atoyan & Aharonian 1999). Relativistic shocks 
should then be formed at the contact interface between 
the ejecta and the wind, providing efficient acceleration 
of electrons. This could result in modulations of the y-ray 
signal with the pulsar spin period (Doppler shifted), which 
would inevitably disappear at times t > Gonf (Eq. 11). 


In Fig. 1^ we show the broad-band spectra of syn¬ 
chrotron and IG radiation predicted by this model at 3 
different times after ejection of the cloud(s) from the ac¬ 
cretion disc: t = 100 s (solid curves), t = Ih (dashed), 
and t = Iday (dot-dashed). For calculations in Fig. 
we have formally assumed that the cloud expands with a 
mean speed ~ 10 ^ cm/s so that the cloud radius at t ~ 1 h 
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Fig. 9. The radiation spectra expected at different times from a compact heavy target with Mc\ = 10^^ g falling 
under a powerful beam of relativistic protons. The fluxes and principal parameters of the cloud, as the radius i?c/, 
the distance Dd from the pulsar, the p — p interaction time tpp, and the cloud temperature Td, are shown. The 
contributions of different radiation processes are shown by solid line (tt^ decays), dashed line (synchrotron radiation), 
dot-dashed line (inverse Compton), 3-dot-dashed line (electron bremsstrahlung), dotted line (’Comptonization’ - the 
multiple Compton scattered radiation in the thermal gas). The heavy solid line corresponds to the total spectrum 
escaping the cloud. Note that the integral fluxes here are integrated ’per decade of the energy interval’, which allows 
to better preserve the information about differential radiation flux as well. 


would still be < cPo/2. The magnetic field behavior in the 
expanding cloud is approximated as B(R) oc with 

a normalisation to Bq = 50 G at Rq = 10^° cm. In such 
high magnetic fields, especially at the very initial stages 
when the cloud is compact, one has to take into account 
that synchrotron losses prevent acceleration of electrons 
effectively beyond the energy Flmax — 50/ ^/bJTG TeV. 
This limitation partially explains an effective suppression 
of TeV radiation in Fig. || at t = 100 s. It is more impor¬ 
tant, however, that the synchrotron self-absorption at the 
initial stages effectively removes from the source all sub¬ 
millimetre and radio photons, which become the principal 
targets for IC production of VHE radiation in such a com¬ 
pact source. 


An important difference between the hadronic and lep- 
tonic compact sources is that the leptonic model does not 
need a high gas density in the source, and the mass load in 
the cloud may be very low. In that case the Coulomb losses 
of relativistic electrons are very low, so the cloud will not 
be heated to very high temperatures when the absorption 
of VHE 7 -rays in the cloud becomes important. Therefore 
only the absorption in the UV radiation field of the 0-star 
may affect the fluxes of VHE y-rays which we do not take 
into account in Fig. For the spectra at t = 100s this 
absorption can be rather strong, depending on the orbital 
phase of the pulsar. Note that at that times the spectral 
modifications are dominated by the absorption of VHE 7 - 
rays in the thermal UV radiation field of the cloud itself. 
At the times t = 1 h, which is the basic timescale for pro- 
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Fig. 10. The spectra of synchrotron (thin lines) and 
IC (heavy lines) radiations from a compact fast cloud 
(‘ejecta’) propagating in the jet(s) of Cen X-3 calculated 
for 3 different times after ejection: t = 100s (solid lines), 
t = Ih (dashed lines), and t = Id (dot-dashed lines). 
The total injection power in shock-accelerated electrons 
Pace = lO^^ergs”^ is assumed. The heavy bar at TeV en¬ 
ergies corresponds to the mean integral flux of TeV y-rays 
detected with the Mark-6 telescope in 1997 (see Chadwick 
et al. 2000) assuming a power-law differential spectrum 
with a = 2.6. 


duction of the pulsed VHE radiation, the source moving 
with Uci > lO'^kms”^ will be at distances far enough from 
the 0-star that the absorption of VHE y-rays is not dra¬ 
matic. Therefore the model predicts rather fast evolution 
of the VHE radiation spectra at the stage of pulsations. 
Note that this prediction is also valid for the hadronic 
compact source model (see Fig. ||). 

Unlike the hadronic model, however, the leptonic jet 
model predicts very signiheant fluxes at later stages of the 
flare evolution, t > 1 d, when the pulsations disappear. It 
is also worth noting that the shape of the spectra in the 
HE y-ray domain in this model are much harder, close to 
the mean spectral flux detected by EGRET in October 
1994 (the hatched region in Fig. 10). Besides that, this 
model can provide an answer to the question why the y-ray 
signals seem to be more persistently detected (although 
at different times) by the Mark-6 telescope in the VHE 
domain than by EGRET at high energies (see dot-dashed 
curve in Fi g. [l(]| ). In this regard note also that the heavy 
bar in Fig. |10| is the mean VHE y-ray flux detected during 
1997, which is higher by a factor of ~ 2 than the mean 
during 1997-1999. 

Finally, in this section we also have to address the ques¬ 
tion of the characteristic frequencies of y-ray pulsations 
that are to be expected. Because the location of the y-ray 
source cannot coincide with the pulsar, one should not 
generally expect that the y-ray pulsations would be co¬ 
incident with the X-ray pulsations. In the framework of 
both models considered in this Section we have to expect 
a difference between those pulsations because of the so 
called ‘double Doppler effect’, which is the frequency shift 


= V — vq due to re-emission/reflection of a pulsed sig¬ 
nal from a target moving with a speed u (see Aharonian 
& Atoyan 19911 ): 




u cos(k'u) — cos(pu) 
c 1 — (u/c) cos(ku) 


( 12 ) 


where u = |u| is the cloud’s velocity relative to the pulsar; 
kTi and pu are the angles that the vector u makes respec¬ 
tively with the y-ray (k) and the relativistic energy out¬ 
flow (p). Then for a compact source moving at high speed 
Uci ^ 10^kms“^ or more favoured by both models above, 
an interpretation of the shift of the Rayleigh power peak 
at 2.400 s from the nominal half-period (second harmonic) 
of the X-ray pulsar Po/2 = 2.410 s would imply that both 
angles ku and pTi were rather close, within < 10°, to 90°. 
In principle, for Gen X-3 where the orbital plane of the 
pulsar is close to the plane of sky this may still be pos¬ 
sible if the jet is produced almost perpendicular to the 
orbital plane. However an interpretation of the y-ray pul¬ 
sations precisely at the X-ray pulsar frequency as reported 


by Vestrand et al. (1997) appears rather problematic. If 
confirmed by future observations, these pulsations could 
be reasonably explained only by a model assuming pro¬ 
duction of HE radiation very close to the pulsar which 
should then address again the problems of the energetics 
and conhnement of such a source. 


4. Conclusions 


The analysis of the data of observations of Gen X-3 accu¬ 
mulated during 23 days of observations of Gen X-3 with 
the University of Durham Mark-6 imaging telescope dur¬ 
ing 1997-1999 conhrms the presence of a y-ray signal at 
the overall signihcance level > 4.5 g in a greement with our 
previous results (Ghadwick et al. 2000 ). The behaviour of 
the signal, which does not drop but rather increases for 
images with enhanced brightness (see Table 1), indicates 
that the spectra of VHE y-rays from Gen X-3 are proba¬ 
bly much harder than the spectra of ambient cosmic rays 
with OCR ~ 2.7. 

No indications for any correlation of the VHE signal 
with the orbital motion of the pulsar, including the pulsar 
eclipse phase, is found. Given that a similar result is found 
by EGRET for HE y-rays detected from Cen X-3, this 
effect excludes any close vicinity of the pulsar as a possible 
site for y-ray production. 

The analysis of the shower arrival times does not re¬ 
veal any statistically significant peak of Rayleigh powers 
in most of the data, except for 1 observation day, the 21 
Feb 1999, showing a strong peak with an estimated prob¬ 
ability of occurence of such a peak by chance, after taking 
into account the large number of IFF trials, of significantly 
below 10“^. The analysis using Bayesian statistics results 
in an even lower final probability, pBay — 7.5 x 10“^. 

It is indicative that the peak powers in both statistics 
are found in the data after application of soft image cuts. 
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whereas they completely disappear when the full (hard) 
image cuts maximizing significance in the overall data of 
observations are applied. This behaviour is just what one 
should expect assuming that the modulated signal in the 
data of 21 Feb 1999 is really connected with VHE y-rays. 
The pulsed signal significantly increases with the thresh¬ 
old energy of the events, which indicates that the spectrum 
of pulsed events is hard as well. 

The position of modulations found in 21 February 1999 
data is blueshifted from the nominal half period of the X- 
ray pulsations by bvjv ~ 3.7 x 10“^. Such shifts should 
be reasonably expected in any model where the y-rays are 
produced far away from the pulsar. Moreover, on the basis 
of the theoretical modelling in this paper, which suggests 
high speeds of compact source(s) needed for production 
of pulsed VHE y-ray emission, we would generally expect 
that y-ray pulsations could show even larger shifts, up to 
\/^v\/vq ~ 0.3 in the case of ejecta moving with the speed 
< 0.3c such as in SS 433. 

The principal models for y-ray production in Cen X- 
3 should be able to produce VHE radiation at distances 
significantly beyond the orbit of the pulsar, effectively at 
i? > 3 X 10^^ cm. Meanwhile, high energies needed for the 
y-ray fluxes in this binary can be provided only by the 
gravitational potential of the neutron star. Confirmation 
of the y-ray fluxes from Cen X-3 by future observations 
would therefore confirm production of powerful jets by the 
inner accretion disc of this X-ray binary. 

Except for the phenomenon of y-ray pulsations (which 
obviously needs confirmation by forthcoming y-ray detec¬ 
tors), an acceptable interpretation of the average fluxes 
currently reported in HE and VHE y-ray regions can in 
principle be provided in the framework of both spatially 
extended and compact source models, which may also 
‘co-exist’ operating together. The principal difference be¬ 
tween these 2 model approaches is that the extended y-ray 
source implies quasi-stationary emission on times scale of 
at least several days, and probably even weeks, whereas 
both hadronic and leptonic compact source models predict 
fast evolution of the y-ray fluxes on time scales of a few 
hours, and with the ‘disappearance’ of the source in a few 
days. Another informative feature could be the behaviour 
of the spectral fluxes at both high and very high energies. 

An important prediction of both compact source mod¬ 
els is that pulsed y-ray emission can only be episodic, with 
a typical duration of no more than a few hours. For y- 
rays of A ~ 100 GeV which are to be produced at large 
distances this is practically a model-independent require¬ 
ment. However, y-rays with E < 10 GeV can escape from 
the binary, even being produced relatively (but not very!) 
close to the X-ray pulsar. At these energies, therefore, de¬ 
tection of pulsed y-rays in the eclipse phase of the X-ray 
pulsar would be very informative. 

We can expect that future observations of Cen X-3 and 
other X-ray binaries with forthcoming sensitive y-ray de¬ 


tectors in the > 50MeV (GLAST) and > 50 GeV (HESS, 
CANGAROO-HI, VERITAS) domains will provide large 
amounts of key information about high energy processes 
in these powerful galactic objects. 
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